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Abstract. A classical theorem of MacMahon states that the number of lozenge tilings of any centrally 
symmetric hexagon drawn on the triangular lattice is given by a beautifully simple product formula. In 
this paper we present a counterpart of this formula, corresponding to the exterior of a concave hexagon 
obtained by turning 120° after drawing each side (MacMahon's hexagon is obtained by turning 60° after 
each step). 



1. Introduction 



MacMahon's classical theorem [16] on the number of plane partitions that fit in a given 
box (see [1] [17] [2] [15] [18] [13] for more recent developments) is equivalent to the fact that 
the number of lozenge"^ tilings of a hexagon of side-lengths a, b, c, a, &, c (in cyclic order) 
on the triangular lattice is equal to 

H(a)H(&)H(c)H(a + 6 + c) 

H(a + 5)H(a + c)H(5 + c)' ^^'^^ 
where the hyperfactorials H(n) are defined by 

H(n) :=0!l!---(n-l)! (1.2) 

(see Figure 1.1 for an example). 

The hexagon is obtained by turning 60° after drawing each side. If instead, after 
drawing each side one turns 120° (the other natural amount of turning on the triangular 
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Figure 1.1. Hexagon with a = 7, Figure 1.2. Shamrock with m = 2, 
6 = 6, c = 8. a = 5, 6 = 7, c = 6. 



lattice), one obtains a shape of the type illustrated in Figure 1.2; we will call such a shape 
a shamrock. 

Our results concern the exterior of a shamrock (as its interior has no lozenge tilings). 
Let S{a,b,c,m) be the shamrock whose central equilateral triangle has side-length m, 
while its top, bottom left and bottom right lobes are equilateral triangles of side- lengths 
a, b and c, respectively; denote its exterior by S*{a, b, c, m). 

We define the ratio of the number of tilings of the exteriors of the shamrocks S{a, b, c, m) 
and S{a-\-b+c, 0, 0, m) as follows. Let H]^{a, b, c, m) he the hexagonal region of side-lengths 
alternating between N+a + b+c and N+a + b+c+rn (the top side being N + a + b+c), and 
having the shamrock S{a, 6, c, to) removed from its center (to be precise, Hi^{a, b, c, m) is 
the region SCN,N,Nici, b, c, to) described in the next section). Then we define 

M(5*(a,&,c,TO)) M{HN{a,b,c,m)) 
M{S*{a + b + c,0,Q,m)) ' n^oo M{HN{a + b + c,0,0,m))' 

where for a lattice region R, M(i?) denotes the number of lozenge tilings of R. 
The dual MacMahon theorem we obtain in this paper is the following. 

Theorem 1.1. For any non-negative integers a, b, c and m we have 

M{S*{a,b, c,to)) _ H(a) H(6) H(c) Rja + b + c + m) H(to)^ 
M{S*{a + 6 + c, 0, 0, to)) ~ H(a + to) H(6 + to) H(c + m)R{a + b + c) 

= P{a,b,m) P{a + b,c,m). (1.4) 



In the special case when m = a + b + c, this can be written simply as 
M{S*{a,b,c,a + b + c)) 



M{S*{a + b + c,0,0,a + b + c)) 



P{a,b,c) P{a + b,b + c,c + a). (1.5) 



Formula (1.5) is illustrated geometrically in Figure 5.1 (see Section 5). Theorem 1.1 
will follow as a consequence of a more general result (see Theorems 2.1 and 2.2), which 
we describe in the next section. The purpose of Section 3 is to establish enumeration 
formulas for so-called "magnet bar regions," which will provide the base cases for the 
inductive proof of Theorems 2.1 and 2.2 given in Section 4. Finally, in Section 5, we show 
how Theorem 1.1 follows from Theorems 2.1 and 2.2. There, we also explain how the 
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Figure 2.1. C6,8,4(2). Figure 2.2. C5,8,4(2). 

results of this paper relate to work of the first author on correlation of holes "in a sea 
of dimers," and we comment on the physical interpretation of our results. We close the 
paper by briefly highlighting the achievements of this paper when compared to the earlier 
paper [3], and what perspectives they offer. 

2. Precise statement of results 

The family of regions whose lozenge tilings we enumerate in this paper is a generalization 
of the following family, introduced in [3]. Consider hexagons of sides ri;, y + m, z, x + m, 
y, z + in (in clockwise order, starting from top) with an equilateral triangle of side m 
removed from its center (see Figures 2.1 and 2.2 for examples). This triangle is called the 
core, and the leftover region, denoted Cx,y,z{m), a cored hexagon. 

To define Cx^y.zi'm) precisely, we need to specify what position of the core is the "central" 
one. Let s be a side of the core, and let u and v be the sides of the hexagon parallel to it. 
The most natural definition would require that the distance between s and u is the same 
as the distance between v and the vertex of the core opposite s, for all three choices of s. 

However, since the sides of the core have to be along lines of the underlying triangular 
lattice, it is easy to see that this can be achieved only if x, y and z have the same parity 
(Figure 2.1 illustrates such a case); in that case, we define this to be the position of the 
core. On the other hand, if for instance x has parity different from that of y and z, the 
triangle satisfying the above requirements would only have one side along a lattice line, 
while each of the remaining two extends midway bcitwccn two consecutive lattice lines (this 
can be seen from Figure 2.2). To resolve this, we translate this central triangle half a unit 
towards the side of the hexagon of length y, in a direction parallel to the side of length x, 
and define this to be the position of the core in this case (see Figure 2.2). 

The family of regions we will be concerned with in this paper is a generalization of cored 
hexagons, corresponding to the case when the core is not just a triangle, but a shamrock. 

Given non-negative integers a, h and c, construct our region as follows. Start with 
the cored hexagon Cx,y,z{m), and "push out" the six lattice lines containing its sides as 
follows: the top, southwestern and southeastern sides a, b and c units, respectively, and 
the bottom, northeastern and northwestern sides b+ c, a + c and a + b units, respectively. 
Enlarge also the core by adding to it downpointing equilateral triangles of sides a, b and 
c that touch the original core at its top, left and right vertex, respectively. The hexagon 
formed by the pushed out edges, with the enlarged, shamrock-shaped core taken out of it, 
is called an S-cored hexagon, and is denoted by SCx,y,z{(i,b,c,m) (see Figures 2.3 and 2.4 
for two examples; the dotted lines indicate the underlying cored hexagons). 

The main result of this paper — from which Theorem 1.1 will follow easily — is the 
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Figure 2.1. 6'C6,8,4(3, 1,2,2). 



Figure 2.2. 5C5,8,4(3, 1,2,2). 



H(n) := 



X 



exact enumeration of the lozenge tilings of S-cored hexagons. This is contained in the 
following two results. In order to formulate them, it is convenient to extend the definition 
of hyperfactorials to half- integers (i.e., odd integers divided by 2): 

YVkZn + l)i for n a positive integer, 
rife^o^ r(fc + i), for n a positive half-integer, 
where F denotes the classical gamma hmction. 

Theorem 2.1. Let x, y, z, a, b, c and m be nonnegative integers. If x, y and z have the 
same parity, we have 

nv./o^ / , NX H(m)'^ H(a) H(6) H(c) 

H(to + a) H(m + b) H(m + c) 
H(2i±H + m + a + 6) H(^ + m + a + c) H(4^ + m + 5 + c) 
H(^ + m + c) H(^ +m + b) H(2i±^ + m + a) 
H(^ + c) H(^ + b) H(^ + g) 
H(^+a + 6) H(2±2+a + c) H(^ + 6 + c) 
H(a; + 771 + a + 6 + c) H(j/ -|-7Ti-|-a-|-6-|-c) 
H(a; -\-y + m + a + b + c) H(a; -|-z-|-m + a-|-6-|-c) 
H(z + m + a + 6 + c) H(a; + y + z + m + a + b + c) 
B.{y + z + m + a + h + c) 
H(r^^1 + m + a + 6 + c)H([^±f±^J + m + a + 6 + c) 
H(^ + m + a + 6 + c) H(^ + m + a + 6 + c) H(^ + m + a + 6 + c) 

H(rf]H(Lfj)H(rf]) 

H([f] +2i+Sd±+£)H([fJ + ™+l^h'!+£)H([|] +Hi+2d±+£) 

H(LfJ)H(rf1)H(LfJ) 

H([f J + "'+"+''+" ) H(rfl + m+a+b+c ^H([fJ + 22+°±b+c) 

H C"'+"+''+'- )-' n(^ + "'+"-''-'• )^ ik;^ + "'+"+''+'• ii(ii±- + 

H([^±|+^] -h rn+a+b+c -)H([£±|±£J + . r»+«+b+c ) H(^) H(^) H(^) ' 



Theorem 2.2. Let x, y, z, a, b, c and m be nonnegative integers. Ifx has parity different 
from the parity of y and z, we have 



M{SCx,y,z{a,b, c,m)) 
H(L^J + m + a4 



H(m)-'^ H(a) H(&) H(c) 
~ H(m + a) H(m + b) H(m + c) 
b) +m + a + c) H(^ 



TO ■ 



b + c) 



TO + c) H([2^J +m + b) H(^+m + a) 
H(r^]+c) H(L^J+6) H(^+a) 
H(L^J +a + 6) H(r^l +a + c) H(^+6 + c) 

H(,T + TO + a + 6 + c) H(2/ + TO + a + 6 + c) 
H(a; + y + m + a + b + c} H(x + 2 + TO + a + 6 + c) 
H(z + TO + a + 6 + c) H(a; + y + z + m + a + b + c) 
Il{y + z + m + a + b + c) 
H(r^±|±^l +m + a + 6 + c)H(L^±|±^J+m + a + 6 + c) 



H([^J +m + a + b + c) H([^l + 

H(rfi)H(Lij)H(rii) 

H([|] + Hi+£±ii+£)H(LfJ + "+°+fe+^ )H(rf 
H(LfJ)H(r|1)H(LfJ) 

H([|J + Hi+^hkhc) H([|] + 12+^^) H([f J + m+a+b+c ^ 



■a + b + c) H(^ + m + a + b + c) 



2 ^ 



1 + 



2 I ' 2 ^-'-'VL 2 J ' 2 J -^-^W 2 1^ 

H([2±2] + r»+a+b+c )H([:£±£J + "'+°+^+° )H(^ + '»+a+fe+c )2 



H([^J)H(M±£) 



(2.2) 



This represents a common generalization of MacMahon's formula (1.1) and Theorem 1.1. 

It shows that, if one regards the hexagon as being the right outer boundary to consider 
on the triangular lattice for the corresponding region to have a number of lozenge tilings 
given by a simple product formula, then a good inner boundary is the shamrock. 

Note also that Theorems 2.1 and 2.2 generalize the main results of [3], by introducing 
three new parameters to the geometry of the core (the sizes of the three lobes of the sham- 
rock). This results in a new, four parameter generalization of MacMahon's theorem (1.1). 



3. Magnet bar regions 



Our proof of Theorems 2.1 and 2.2 will use the exact enumeration of the following 
regions. 

Let X, y, a, b, c and m be nonnegative integers. Consider the hexagonal region with two 
removed equilateral triangles illustrated by Figure 3.1, where the top side lengths of the 
hexagon are, clockwise from top, x + c, y + m, a + b + c, x + m,y + c, a + b + m, and the sides 
of the two removed triangles are m (for the triangle touching the northwestern hexagon 
side; note that the lengths of the portions of that side above and below this triangle are 
a and b, respectively) and c (for the central triangle). We call such a region a magnet bar 
region and denote it by Bx^y{a, b, c, to). 

The main result of this section is the following. 
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x+c 




Theorem 3.1. For any nonnegative integers x, y, a, b, c and m, we have 
H(m)2 H(a) H(6) H(c) H(m + a + h + c) 



M{Bx^y{a,b,c,m)) = 



H(m + a) H(m + b) H(m + c) 

ii{x + m + a + c) H(y + m + b + c) R{x + y + c) 

Rix + y + m + c) H(x + a + c) H(j/ + 6 + c) 

H(a; + y + m + a + 5 + c) H(a:) H(2/) 



H(a; + m + o + 6 + c) H(j/ + m + a + 6 + c) H(3; + y) 



(3.1) 



The proof of the above result (as well as the proof of Theorems 2.1 and 2.2) is based on 
Kuo's powerful graphical condensation method (see [14]). For ease of reference, we state 
below the particular instance of Kuo's general results that we need for our proofs (which 
is Theorem 2.1 in [14]). 

Theorem 3.2 (Kuo). Let G = (Vi, V2,E) be a plane bipartite graph in which \Vi\ = IV2I. 
Let vertices a, /3, 7 and S appear cyclically on a face of G. If a,^ € Vi and P,d € V2, 
then 



M(G) M(G - {a, /3, 7.5}) = M(G - {a, /3}) M(G - {7, S}) + M(G - {a, S}) M(G - 7}). 

(3.2) 

Proof of Theorem 3.1. We prove (3.1) by induction on x + y + b. Our base cases will 
be the situations when a; = 0, y = or 6 = 0. If 5 = 0, note that the lozenge tiling is 
forced in an m x (y + c) rectangle along the southwestern edge of the magnet bar region 
(Figure 3.2 illustrates this for x = 3, y = 1, a = 4, c = 3, m = 2). Upon removing these 
forced lozenges, the leftover region is one whose tilings are enumerated by a formiila due to 
Cohn, Larsen and Propp (see Proposition 2.1 in [10]). Indeed, regard the leftover region 
— shown in bold contour in Figure 3.2 — as being obtained from a hexagon by cutting 
out from it a triangle of side c resting on its southwestern edge. It is readily seen that this 
region has the same number of lozenge tilings as the region obtained from the hexagon 
by cutting out c consecutive unit triangles resting on the southwestern edge (this is due 
to forced lozenges in any tiling of the latter). However, regions of this type have their 
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x+m 



m 



Figure 3.2. The case 6 = 0. Figure 3.3. The case x = 0. 

number of lozenge tihngs given by [10, Proposition 2.1]. It is routine to check that the 
resuhing formula agrees with the 6 = specialization of the expression on the right-hand 
side of (3.1). 

If X = 0, it turns out that tiling the magnet bar region with lozenges is equivalent to 
tiling two disjoint, hexagon shaped subregions (this is illustrated in Figure 3.3; the hexag- 
onal subregions are determined by the two dashed lines; that one of them is interrupted 
by a shaded lozenge should be ignored at this point). 

To see this, consider for instance the top hexagon Hi, obtained by extending upward 
the southeastern side of the triangle of side c (shown in Figure 3.3 by a dashed line). Let 
T be a lozenge tiling of the magnet bar region. We claim that there is no lozenge in T 
that straddles the southeastern side of Hi. Indeed, suppose towards a contradiction that 
T contained such a lozenge L (pictured in gray in Figure 3.3). Following upward from L 
the lozenge in T above L, then the lozenge in T above that, and so on, one obtains a path 
of lozenges that eventually must end on the upper side of the magnet bar region. However, 
the top side of the cut out triangle of side c generates in a similar fashion c more paths 
of rhombi, which must also end at the top side of the magnet bar region. Since the latter 
side has length c and the above c + 1 paths of rhombi are necessarily nonintersecting (as 
they are all part of the tiling T) , this provides the contradiction that proves our claim. 

Thus the top dashed line cannot be crossed by any lozenge. It follows that the rhombus 
cut out from the eastern corner by the two dashed lines is forced to be part of each tiling, 
and can therefore be removed without affecting the number of tilings of the magnet bar 
region. The leftover region is the union of two disjoint hexagons. Using (1.1) it follows 
that for a; = we have 

M(Bo,3/(a, b, c, m)) = P{a, c, m) P{b, y + c,m). (3.3) 

It is apparent from (1.1) that the right-hand side above agrees with the x = specialization 
of the right-hand side of (3.1). The case y = follows similarly. This completes the 
verification of the base cases of our induction. 

The induction step is based on a convenient application of Kuo's graphical condensation 
stated in Theorem 3.2. 

Suppose x,y,b > 1, and assume that (3.1) holds for all values of the parameters for 
which the sum of the x-parameter, ^/-parameter and 6-parameter is less than x + y + b. 
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Figure 3.4. The recurrence for the regions Bx,y{a,b,c,m). 



Let G be the planar dual graph** of the region Bx^y{a, b, c, m). Choose the vertices a, /3, 
7 and 6 as indicated in Figure 3.4, where a is the leftmost black unit triangle^, and j 
and 5 arc the next black unit triangles as one moves counterclockwise from a (Figure 3.4 
corresponds to the case x = 3, y = 1, a = 4, 6 = 1, c = 3, m = 2.) Then (3.2) states 
that the product of the number of lozenge tilings of the two regions on top is equal to the 
product of the number of lozenge tilings of the two regions in the middle, plus the product 
of the number of lozenge tilings of the two regions on the bottom. After removing the 
lozenges forced by the unit triangles a, /3, 7 and 6 indicated in Figure 3.4, the leftover 



*By the planar dual graph of a region on the triangular lattice we understand the graph whose vertices 
are the unit triangles inside the region, and whose edges connect vertices corresponding to unit triangles 
that share an edge. 

®By this we understand, of course, the vertex of the planar dual graph corresponding to that unit 
triangle. 
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regions are, in all six instances, magnet bar regions. More precisely, we obtain 

M{Bx,y{a, b, c, to)) M(i?a;_i,y (a, 6- 1, c, m)) = M{Bx^y{a,b- 1, c, to)) M{B^_i^y{a, b, c, to)) 

+ M{B^_i^y+i{a,b-l,c,m))M{B^^y_i{a,b,c,m)). (3.4) 

Note that all magnet bar regions in the above equation except the first one have the sum 
of their x-, y- and 6-parameters strictly less than x + y + b. By the induction hypothesis, 
these five magnet bar regions have their number of lozenge tilings given by (3.1). It is 
readily checked that substituting these formulas into the above equation one obtains, after 
simplifications, that M{Bx,y{a,b,c,m)) equals precisely the expression on the right-hand 
side of (3.1). This completes the induction step, and hence the proof. □ 

4. Proof of Theorems 2.1 and 2.2 

As we have already mentioned, our proof of Theorems 2.1 and 2.2 is also based on 

Kuo's graphical condensation method. The argument we use is an induction that proves 
Theorems 2.1 and 2.2 simultaneously. The previous section's result on magnet bar regions 
is needed for the base cases of our induction. 

In our proof it will be essential to be familiar with various distances one can naturally 
consider within a given ^-cored hexagon. These are displayed in Figures 4.1 and 4.2 (note 
that the unit used is the distance between two consecutive lattice lines). 

Proof of Theorem,s 2.1 and 2.2. We prove formulas (2.1) and (2.2) by induction on 
X + y + z. The base cases are the instances when x = y = Q ov z = 

Consider first the case when x, y and z have the same parity. Due to symmetry, 
it is enough to verify the instance when z = 0. In that case, the 5-corcd hexagon 
SCx,y,o{a, b, c, to) looks as illustrated in Figure 4.3. Then the hexagon H cut out from the 
S-coied hexagon by extending sides of the a- and 6-lobes of the shamrock core as indicated 
in Figure 4.3 (see the dotted lines in the figure) must be internally tiled in each tiling of 
SCx,y,Q{a, b, c, to). Indeed, this follows by the same argument we used to show that the top 
hexagonal subregion in Figure 3.3 is internally tiled, except we need to apply it now for 
both dashed-linc cuts (the argument applies because the northwestern and southeastern 
sides of H have the same length, to). 

Since H is tiled internally, it follows that the rhombus that fits in the top right corner 
of SCx.y.o{a,b, c,m) and rests on both H and the a-Iobc of the shamrock is forced to be 
tiled as shown. The same is true for the analogous rhombus fitting in the bottom left 
corner. Note that the region obtained from SCx,yfl{a,b,c,m) by removing H and these 
two rhombi is precisely the magnet bar region B:^ :K{a,b,c,m) (see Figure 4.1; since we 
are in the case when x, y and z have the same parity and we are assuming z = 0, both x 
and y are even). It follows that we have 

M{SCx,y,oia, b, c, to)) = M{H) M{B.k (a, b, c, to)) 

= P (to, I + 6, 1 + a) M{B.^y_ (a, 6, c, to)) (4.1) 

(for the second equality we used Figure 4.1 to read off the side- lengths of H). Substituting 
the values of P(to, 1 + 6, 1 + a) and M(i?|^£ (a, b, c,m)) given by formulas (1.1) and (3.1) in 
the right-hand side above, it is apparent that it becomes precisely the z = specialization 
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x+a+b+c 




x+m 



Figure 4.1. Distances within the 5-cored hexagon SCx,y,z{a, b, c, m) 
when X, y and z have the same parity. 



x+a+h+c 




x+m 



Figure 4.2. Distances within an S'-cored hexagon SCx,y,z{o-,b,c,m) 
when X has opposite parity to y and z. 

of the expression on the right-hand side of (2.1). This checks the base case of the induction 
when X, y and z have the same parity. 

When X, y and z have mixed parities, we can assume without loss of generahty that x 
has parity opposite to the parities of y and z. There are now two inequivalent base cases 
to check, namely .t = and z = (y = is equivalent to the latter). Consider first the 
case when x = (note that our assumptions imply then that y and z are odd). Then the 
S'-cored hexagon SCo^y^z{0',b,c,m) looks as illustrated in Figure 4.4. 

Note that, by the argument we used in the paragraph before (4.1), the hexagonal 
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Figure 4.3. The S'-cored hexagon 5(70,4,4(3, 1,2,2). 




Figure 4.4. The region S'Co,5,3(3, 1,2, 2). Figure 4.5. The region 5(75,4,0(3,1,2,2). 

subregion determined by the dashed Unes in Figure 4.4 must be internally matched in any 
lozenge tiling of SCo^y^z{a,b,c,m). This implies that the two rhombic regions at the left 
and right corners of the S-coied hexagon whose lozenge tilings are indicated in Figure 4.4 
are tiled as shown in each tiling of 5(7o,j/,z(a, b, c, m). The region left over from the S-coved 
hexagon after removing these two rhombic regions is the disjoint union of the hexagonal 
subregion mentioned earlier in this paragraph, and a magnet bar region. It follows from 
the distances indicated in Figure 4.2 that the three values for the hexagon side-lengths are 
m, + c and + b, while the magnet bar region is By^ z±i {b, c, a, m). Thus we have 

M(5Co,j,,z(a, 6, c, to)) = P (^TO, + c, + b^ M{By-i .±i{b, c, a, to)). (4.2) 

Using formulas (1.1) and (3.1) one readily sees that the above equation agrees with the 
a; = specialization of (2.2). 

The remaining base case to check is z = 0, x odd and y even. This follows by the same 
arguments as the other two base cases. Figure 4.5 shows how M{SCx,yfl{a, b, c, m)) is the 
product of the number of tilings of a hexagon and a magnet bar region. Their precise 
dimensions can be deduced from Figure 4.2. One obtains 

M{SC:c,y,oi(i^ b, c, to)) = P (m, + &, I + M(Bx+i _^ (a, b, c, to)). (4.3) 
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Figure 4.6. The recurrence for the regions SCx,y,z{0,0,0,Tn), x of opposite parity to y and 

Using (1.1) and (3.1), this is again readily seen to agree with the z = Q speciahzation of 
(2.2). This conchides the verification of all the base cases we need for our induction. 

The induction step is based on Kuo's graphical condensation (stated in Theorem 3.2). 
More precisely, we use graphical condensation to obtain recurrences for the number of 
lozenge tilings of the S'-corcd hexagons SCx.y.z{o,,b,c,m), valid for x,y,z > 1. and then 
verify that the right-hand side of the claimed formulas (2.1) and (2.2) satisfy the same 
recurrences. 

We will need two different recurrences, one for the case whcin ,t, y and z have the 
same parity, and one for the mixed parity case. Both will follow by applying graphical 
condensation in the same fashion as in the proof of Theorem 3.1, namely by choosing the 
four removed unit triangles to be along two sides of the outer boundary of the S'-corcd 
hexagons, in the pattern shown in Figure 3.4. In particular, only the outer boundary of 
the 5-cored hexagons will change in the resulting subregions, while the shamrock core 
remains intact. Due to this, it suffices to discuss how we obtain our recurrences in the 

12 



case when a = b = c = 0. 

Consider first the case when x, y and z have mixed parities (an instance of this is 
illustrated in Figure 4.6). Without loss of generality we may assume that the parity of x is 
different from that of y and z. When applying Kuo's graphical condensation as described 
in the previous paragraph, it is of crucial importance in which corner of the outer hexagon 
one places the pattern of the four removed unit triangles. Since in SCx,y,z{0,0,0,m) 
(which is the same as the corcHl hexagon C'x.y.z{m) of [3]) the core is just to the left of the 
true central position, the correct choice will be to place the pattern of the four removed 
unit triangles in the right corner of the outer hexagon (sec Figure 4.6). 

Let G be the planar dual graph of the region SCx,y,z{0, 0, 0, m). Choose the vertices a, 
/3, 7 and 5 as indicated in Figure 4.6, where a is the bottom black unit triangle, and /3, 7 
and 6 are the next black unit triangles as one moves upwards (Figure 4.6 corresponds to 
the case .x = 4, y = 7, z = 3). Then (3.2) states that the product of the number of lozenge 
tilings of the two regions on top is equal to the product of the number of lozenge tilings 
of the two regions in the middle, plus the product of the number of lozenge tilings of the 
two regions on the bottom. After removing the lozenges forced by the black unit triangles 
indicated in Figure 4.6, the leftover regions are, in all six instances, S'-cored hexagons. 
(Note that this would not be the case if we applied graphical condensation with the 
pattern of rcmovcHl triangles placed in a different corner — some of the resulting positions 
of the core would not be central!) The precise parameters of these S'-cored hexagons can 
be visually extracted from Figure 4.6 (this is easier to do in the case o = 6 = c = 0; this 
is why we have reduced the proof of the recurrences to this case). 

Indeed, the leftover region on the top right of Figure 4.6 is SCx,y-i,z-i{0, 0, 0, m). The 
central left region, after clockwise rotation by 120° followed by a reflection across the verti- 
cal, becomes SCy^x,z-i{0, 0, 0, m). The central right region, after counterclockwise rotation 
by 120° followed by a reflection across the vertical, becomes SCz,y-i,x{Q, 0, 0, m). Similarly, 
the bottom two leftover regions are SCx-i,y,z{Q, 0, 0, m) and SCx+i,y-i,z-i{^, 0, 0, m), re- 
spectively. We obtain 

MiSCx,y,z{0, 0, 0, m)) MiSCx,y-i,z-i{0, 0, 0, m)) 

= M{SCy,x,z-i{0, 0, 0, m)) M{SCz,y-i,x{Q, 0- 0, m)) 

+ M{SCx-i,y,ziO, 0, 0, m)) M{SCx+i,y-i,z-iiO, 0, 0, m)). (4.4) 

As explained above, this yields a recurrence for any nonnegative values for the sizes a, b 
and c of the lobes of the shamrock core. Taking into account the two rotations followed 
by reflections that we needed to consider when converting Figure 4.6 into recurrence (4.4), 
one sees that the resulting recurrence for arbitrary a, 6, c is 

M{SCx,y,z{0', b, c, m)) M{SCx,y-i,z-i{o-^ b, c, m)) 

= M{SC'y^x,z-i{b, a, c, m)) M{SCz,y-i,xic, b, a, m)) 

+ M{SCx-i,y,z{a, b, c, to)) M{SCx+i,y-i,z-i{a, b, c, to)), 

x,y,z> 1, X of opposite parity to y and z. (4.5) 
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Figure 4.7. The recurrence for the regions SCx,y,z{0,0,0,'m), x, y and z of same parity. 

In the remaining case when x, y and z have the same parity, there is no "singled 
out" corner of the outer hexagon in which to place the pattern of removed unit trian- 
gles when applying condensation any of the six choices gives a rcciirrcncc for S'-cored 
hexagons. In order to be closer to the previous case, we choose again to fit this pat- 
tern in the rightmost corner; see Figure 4.7, which illustrates the case x = 5, y = 7, 
z = 3. Then (3.2) states that the product of the number of lozenge tilings of the two 
regions on top is equal to the product of the number of lozenge tilings of the two regions 
in the middle, plus the product of the number of lozenge tilings of the two regions on 
the bottom. After removing the lozenges forced by the black unit triangles indicated in 
Figure 4.7, all the leftover regions are ^-cored hexagons. Their precise parameters can 
be easily deduced from Figure 4.7. Indeed, the region on the top right, after reflection in 
the vertical, becomes SCx^z-i.y-ii(>, 0, 0, m). After counterclockwise rotation by 120°, the 
left central region becomes SCz-i,y,x{0,0,0,m). The right central region is transformed 
into SCy-i^z^x{0, 0, 0, m) by clockwise rotation by 120°. The bottom left region, when re- 
flected across the vertical, becomes SCa,-! ^ ,^(0, 0, 0, m). Finally, the bottom right region 
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is precisely SCx+i,y-i,z-i{0, 0,0,1^1). We obtain 
M(5a,^,,(0, 0, 0, m)) M(5C,,,_i,y_i(0, 0, 0, m)) 

= M(5C,_i,j,,,(0, 0, 0, m)) M(5Cj,_i,,,,(0, 0, 0, m)) 

+ M(5C,_i,,,^(0, 0, 0, m)) M(5C,+i,^_i,,_i(0, 0, 0, m)). (4.6) 

By the above diseussion, taking into aecount the symmetries we needed to apply to the 
regions in Figure 4.7 in the previous paragraph, this extends to arbitrary values of a, b 
and c. The resulting recurrence is 

M{SC^,y^,{a, b, c, m)) M(S'a,^_i,j,_i(a, c, b, m)) 

= M(S'C^_i,a,,j,(c, a, b, to)) M(S'Cj,_i,^,^(6, c, a, to)) 

+ M{SCx-i,z,v{a, c, b, to)) M{SCx+i,y-i,z-i{a, b, c, to)), 

x,y, z > 1, a;, y and z of same parity. (4.7) 

Let x,y,z > 1 and a, b, c and to be fixed nonnegative integers, and assume that (2.1) 
and (2.2) hold for all values of the parameters for which the sum of the a;-parameter, 
y-parameter and z-parameter is less than x + y + z. We need to show that (2.1) and (2.2) 
hold also for x, y, z, a, b, c and to. 

Note that recurrences (4.5) and (4.7) express M{SCx,y,z{a,b,c,rn)) in terms of five 
quantities of the form M{SCxi ^yi ^z' (o-, b, c, m)) , with x' + y' + z' < x + y + z for each of 
them. By the induction hypothesis, formulas (2.1) and (2.2) hold for each of these five 
quantities. To complete the proof, it suffices to show that the expression obtained for 
M{SCx,y,z{o-,b,c,'m)) when substituting the formulas given by (2.1) and (2.2) for these 
five quantities is equal to the right-hand side of (2.1) (for recurrence (4.5)), respectively 
the right-hand side of (2.2) (for recurrence (4.7)). Defining 



Rx,y,z{0;b,c,m) := < 



right-hand side of (2.1), x, y, z of same parity, 

(4.8) 

right-hand side of (2.2), x of parity opposite to y and z, 
this is equivalent to verifying that 

Rx,y,zia, b, c, to) Rx,y-i,z-i{a, b, c, to) = Ry^x,z-iib, a, c, to) Rz^y-i,xic, b, a, to) 

+ -Rx-i,!/,z(o, b, c, to) Rx+i,y-i,z-i{a, b, c, to), 
x,y,z> 1, X of opposite parity to y and z (4.9) 

and 

Rx,y,z{(^^ ^ c, to) i?x,2-i,j/-i(a, c, b, to) = Rz-i,x,y{c, a, b, to) Ry-i^zA^^ c, a, to) 

+ Rx-i,z,y{a, c, b, to) Rx+i,y-i,z-i{a, b, c, to), 

x,y,z> 1, X, y and z of same parity. (4-10) 
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Wc verify (4.10) first, so assume that x, y and z have the same parity. We compute the 
quotient of the right-hand side by the left-hand side in (4.10). Using the "cancellation 
rule" 

H(r.T + yl) H(Lx + yJ) 



H(rx-il+y)H(Lx-iJ+2/) 



^T{\x + y-]), (4.11) 



where x is an integer or a half-integer, this quotient can be significantly simplified. Namely, 
after applying it numerous times, it becomes 

T{y + z + a + b-\-c + m)T{x + y + z + a + h + c + m-l) 
V{y + z + a + b + c + m-l) T{x + y + z + a + h + c + m) 

r ([^±^1 ru + h + c + rn) r ( \^] + '-^i^) 

r([^±f±^] +a + b + C + m) r([ W^-l ^ ^ a+fe+c+m -) 

T{x + a-\-b + c + m+l)T{x + y + z + a + b + c + m-l) 



+ 



T{x + a + b + c + m)T{x + y + z + a + b + c + m) 

r(r^l)r([f]+i^±^) 



■r([fl)r(m + ^^±^: 



+a + b + c + m) r([^±f±^] + °+fc+°+'» ) 
+ a + 6 + c + m) r([^^±^±^] + 

Let first x, y, z all be even. Then the expression in (4.12) reduces to 

(y + z + a + b + c + m— 1) 
(x + y + z + a + b + c + m-l) 

{x + a + b + c + m) x/2 



(4.12) 



(x + y + z + a + b + c + m-l) {x + a + b + c + m)/2' 

which indeed equals 1. 

Now we assume that all of x, y, z are odd. In that case, the expression in (4.12) reduces 

to 

(y + z + a + b + c + m-1) ( W^-i ^ a+b+c+m ^ 



2 ' 2 



(x + y + z + a + b + c + m-l) [2=±b±£^ + a + b + c + m) 



^ ix + a + b+C+m) ( .+^+.-1 ^ a+b+c+m ^ 



(x + y + z + a + b + c + m-l) 



- m 



which also equals 1. 



We next turn to verifying (4.9), so assume that x has parity different from the parity 
of y and z. Here, we compute the quotient of the right-hand side by the left-hand side in 
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(4.9). After applying the cancellation rule (4.11) numerous times, this quotient becomes 

T{y + z + a + b + c + m) F{x + y + z + a + b + c + m— 1) 
T{y + z + a + b + c + m-l)T{x + y + z + a + b + c + m) 

r([^±a±^] +a + b + c + m) r([^^±f±i] + "+fe+^+'» ) 



(i-\-b-\-c-\-m ') 
2 ) 



+ 



r([^^l +a + b + c + m) r([ W"-^ ] + 

r(a; + a + b + c + m + l)r{x + y + z + a + b + c + m - 1) 
T{x + a + b + c + to) r(a; + y + z + a + b + c + m) 

r(m)r([fl+i^±^) 



r([fl)r(m + 



^, r(r^±^1 +a + 6 + c + m) r(rm] + .±6±^) 
r(r^^1 +a + b + c + m) r(r^+^] + "+''+-+'" )• ^ ^ ^ 

Note that this happens to be identical to the expression (4.12). Since y + z is even in our 
case, regardless of the particular choice of y and z, and the expression in (4.13) depends 
only on x and y + z (and not on y and z separately!), we may use again the earlier 
computations to deduce that this always simplifies to 1. 

This completes the verification of (4.9) and (4.10), and hence the proof. □ 

Remark 1. The special case a = b = c = Oo{ Theorems 2.1 and 2.2 was the main result 
of our earlier paper [3] (sec Theorems 1 and 2 there). The proofs given in [3] used the 
method of factor exhaustion, and required considerable calculations (and space!). The 
proof we presented above represents in particular a great simplification over the proofs 
in [3]. 
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Remark 2. At the other extreme;, the; speeiahzation x = y — z — oi Theorems 2.1 
and 2.2 provides the following interesting picture (see Figure 4.8). In this case, the S- 
cored hexagon S'Co,o,o(a, b, c, m) is the disjoint union of three hexagons. It follows that 

M(S'Co,o,o(a, c, m)) = P(a, b, m) P{a, c, m) P{b, c, m). (4-14) 

This offers a different way of viewing Theorems 2.1 and 2.2 as a generalization of MacMa- 
hon's theorem (1.1). Namely, they can be regarded as simultaneously generalizing three 
applications of MacMahon's formula. 

5. Deducing Theorem 1.1. Geometric and physical interpretation 

Proof of Theorem 1.1. By (1.3), the explicit definition of the expression on the left-hand 
side of (1.4) is 

M{S*{a,b,c,m)) _ ^.^ M(5'Cjv,Af,jv(a, 6,c, m)) 



M(S'*(a + 6 + c,0,0,m)) ' jv-j-oo M(5Cjv,jv,Jv(a + & + c, 0, 0,m)) " 

In the S-cored hexagons on the right-hand side above, the x-, y- and z-parameters have the 
same parity (being all equal to N). Thus the number of their lozenge tilings is given by the 
formula provided in Theorem 2.1. Applying this formula one obtains, after simplifications, 
that 

H(m)^ H(m + a + & + c) 
M{SCN,N,N{a., 6, c, m)) _ H(m + a) H(m + 6) H(m + c) 



M{SCN,N,N{a + b + c,Q,0,m)) H(a + b + c) 



H(a) H(6) H(c) 

H(7V + m) H(7V + m + a + b) H(A^ + m + a + c) H(7V + m + b + c) 
H(jV + m + g) H(jV + m + b) H(iV + m + c) H(jV + m + a + 6 + c) 

^ H(iV) njN + a + h) H(iV + a + c),H(iV + bT^ 2) 

H(7V + a) H(Ar + fe) H(Ar + c) H(7V + a + 6 + c) 

Recall that, by the Glaisher-Kinkelin formula (see [12]) which gives the asymptotics of 
the Barnes G-function, we have 

0!1! ••• (n- 1)! 

lim ^ , 1. = (5.3) 

n'^-ii (27r)t e-^ ^ 

where ^4 = 1.28242712... is a constant (called the Glaisher-Kinkelin constant). This readily 

implies that 

H(jV) H(iV + a + 5) H(iV + a + c),H(iV + 6 + c) ^ 
jv^ H(7V + a) H(7V + 6) H(7V + c) H(7V + a + 6 + c) ' ^ ' 

As m is fixed, (5.4) shows that both the numerator and the denominator of the second 
fraction on the right-hand side of (5.2) approach 1 as iV — ^ 00. It follows then from (5.2) 
that 

II(//( + a + + c) 

J. M{SCN,N,N{a':b,c,m)) __ H(to -|- a) H(m + 6) H(m -|- c) , , 

M(S'Cjv,Ar,jv(a + 6 + c, 0, 0, m)) ^ H(a + 6 + c)~ " ^ ' 

H(a) H(6) H(c) 

The first equality in (1.4) follows from (5.1) and (5.5). The second equality in (1.4), as 
well as (1.5), follow then using formula (1.1). □ 
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Figure 5.1. Geometric interpretation of Theorem 1.1 
when m = a + b + c; here a = 3, 6 = 1, c = 2. 



The hexagons occurring in the symmetric form (1.5) of Theorem 1.1 (corresponding to 
the case when m = a + b + c) have the foUowing geometric interpretation. Reflect each 
lobe of the shamrock across the vertex of the central triangle touching that lobe. This cuts 
out from the central triangle of the shamrock a hexagonal region, whose side-lengths are, 
clockwise from top, a, b, c, a, b, c. We call this the central hexagon of the shamrock. For 
a = 3,6=l,c = 2, this is indicated by the solid line contour in the center of Figure 5.1. 

Another natural hexagon is obtained as follows. Consider the six pairs of next-to-nearest 
edges of the central hexagon (we call two edges of a hexagon next-to-nearest if they are not 
incident to one another, but arc both incident to another edge of the hexagon). Three of 
these pairs cross at the vertices of the central triangle of the shamrock. Extend the other 
three pairs as well to obtain three more crossing points. Define the medial hexagon of the 
shamrock to be the smallest lattice hexagon containing the above six crossing points. The 
solid line contour on the right in Figure 5.1 illustrates this in the case a = 3, 6 = 1, c = 2. 

Then formula (1.5) can be expressed geometrically as shown in Figure 5.1: The ratio 
of the number of tilings of the exteriors of the two shamrocks is equal to the product of 
the number of tilings of the central and medial hexagons of the shamrock. 

Remark 3. In the language of the first author's earlier series of papers on the correlation 
of holes in a sea of dimers and their connection to two dimensional electrostatics (see 
[4] [5] [6] [7] [8]), the right-hand side of (1.3) is the ratio between the correlation u! defined 
in [4] of the two shamrocks. Therefore (1.4) can be restated as 

^(5(a,6,c, m)) _ H(o) H(6) H(c) Rja + b + c + m) H(m)^ 
uj{S{a -I- 6 -f- c, 0, 0, m)) ~ H(a + m) H(fe -|- m) H(c -|- m) H(a + b + c) 

= P{a,b,m) P{a + b,c,m), (5.6) 
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and (1.5) as 

Lu(S(a,b,c,a + b + c)) „, ^ ^ , ^ 

^ \' ' ' ^ -^^—-=P{a,b,c)P{a + b,b + c,c + a). 5.7 

iu{S{a + b + c,0,0,a + b + c)) v > > y v ; \ j 

In particular, it follows that for all values of the parameters a, b, c and m, the ratio of the 
two correlations on the left-hand side of (5.6) is an integer. There does not seem to be 
any simple a priori reason why this should be so. 

Theorem 1.1 also has an interpretation in the light of the connection to electrostatics 
cited above. According to this, holes on the triangular lattice in a sea of lozenges interact 
(when their correlation is taken, or measured) precisely like charges in two dimensional 
electrostatics. Namely, the correlation of the holes is equal, in the limit of large mutual 
separations between the holes, to the exponential of the negative of the electrostatic energy 
of the two dimensional system of electrical charges obtained by regarding each hole as an 
electrical charge, of magnitude equal to the number of up-pointing unit triangles in the 
hole minus the number of down-pointing unit triangles in the hole. 

In the set-up of Theorem 1.1, we have four triangular holes making up the shamrock 
whose correlation is considered. One readily sees that in the above described connection to 
electrostatics, the lobes correspond to charges of magnitudes a, b and c, while the central 
triangle corresponds to a charge of magnitude — m. The three lobes should then repel 
one another. And indeed, by (5.6) the least likely configuration (i.e., the configuration of 
maximum energy), when one fixes m and the value of the sum a + b + c, is the situation 
when one lobe contains the entire positive charge, and the other two lobes vanish. 

Note that formulas (5.6) and (5.7) are relative correlation results, expressing the ratio of 
the correlation of two shamrocks. We end this remark by explaining how the correlations 
on the left-hand side of (5.7) can be computed separately. 

To do this, since we have formula (5.7), it suffices to determine Lo{S{a -|- 5 -|- c, 0, 0, m)). 
After a straightforward (if somewhat lengthy) calculation, using Stirling's formula and the 
Glaisher-Kinkelin formula (5.3) to analyze the asymptotics of the resulting expressions of 
type (1.1) and (2.1), one obtains that 

/r./ , XX M(SC^^Jm,0,0,m)) Vs"' H(m)4 

w(5(a-|-6-Fc,0,0,m)) := lim — — ^ ' ' ' ' — '-^ ^ , ^^)^ \ (5.8) 

^ ^ ' ' ' " x^cc.M{H{x + m,x + m,x + m)) (27r)™ H(2m) ^ ' 

(here H{x + m,x + m,x + m) denotes the regular hexagon of side x + m). 

Remark 4. We end this section with an interpretation of Theorems 2.1 and 2.2 in terms of 

surfaces analogous to plane partitions, but having a certain defect in the middle. Wc con- 
fine our discussion to the case m = a + b + c (the general case also admits an interpretation 
along these lines, but it is less elegant). 

Recall that a lozenge tiling of a hexagon naturally lifts to a three dimensional surface 
(see e.g. [11]). This lifting works also in the presence of a hole of charge zero (see the 
detailed discussion in [5]). For instance, the tiling shown on the left in Figure 6.1 lifts to 
the surface shown on the right in the same figure (note that this surface has an overhang 
and some precipices in the middle). From this point of view, the case m = a + b -\- c 
of Theorems 2.1 and 2.2 states that the number of such surfaces is given by the explicit 
product formulas (2.1) and (2.2). 

6. Concluding remarks 



We have seen in this paper the great power of Kuo's graphical condensation method, 
when applied to situations in which explicit conjectured formulas can be found (for another 
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Figure 6.1. A tiling of 56,6,6(1, 2, 1,4) and its lifting. 



application of the same method, see [9]). The proofs presented here represent a great 
simplification of the proofs of Theorems 1 and 2 of [3] , which are the special case a = b = 
c = of Theorems 2.1 and 2.2 of this paper. Furthermore, the same method can be used 
to tackle Conjectures 1 and 2 of [3], which seem otherwise quite redoubtable, and also to 
extend them to 5-cored hexagons (details will appear in a separate paper). 
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